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ABSTRACT 

We present an analysis of the X-ray emission of the rapidly rotating giant star YY Mensae observed by Chan- 
dra HETGS and XMM-Newton. The high-resolution spectra display numerous emission lines of highly ionized 
species; Fe XVII to Fe XXV lines are detected, together with H-like and He-like transitions of lower Z elements. 
Although no obvious flare was detected, the X-ray luminosity changed by a factor of two between the XMM- 
Newton and Chandra observations taken 4 months apart (from logLx ~ 32.2 to 32.5 erg s~', respectively). The 
coronal abundances and the emission measure distribution have been derived from three different methods us- 
ing optically thin collisional ionization equilibrium models, which is justified by the absence of opacity effects 
in YY Men as measured from line ratios of Fe XVII transitions. The abundances show a distinct pattern as a 
function of the first ionization potential (FIP) suggestive of an inverse FIP effect as seen in several active RS 
CVn binaries. The low-FIP elements (< 10 eV) are depleted relative to the high-FIP elements; when compared 
to its photospheric abundance, the coronal Fe abundance also appears depleted. We find a high N abundance in 
YY Men's corona which we interpret as a signature of material processed in the CNO cycle and dredged-up in 
the giant phase. The corona is dominated by a very high temperature (20-40 MK) plasma, which places YY 
Men among the magnetically active stars with the hottest coronae. Lower temperature plasma also coexists, 
albeit with much lower emission measure. Line broadening is reported in some lines, with a particularly strong 
significance in Ne X Lya. We interpret such a broadening as Doppler thermal broadening, although rotational 
broadening due to X-ray emitting material high above the surface could be present as well. We use two dif- 
ferent formalisms to discuss the shape of the emission measure distribution. The first one infers the properties 
of coronal loops, whereas the second formalism uses flares as a statistical ensemble. We find that most of the 
loops in the corona of YY Men have their maximum temperature equal to or slightly larger than about 30 MK. 
We also find that small flares could contribute significantly to the coronal heating in YY Men. Although there is 
no evidence of flare variability in the X-ray light curves, we argue that YY Men's distance and X-ray brightness 
does not allow us to detect flares with peak luminosities Lx < 10^' erg s"' with current detectors. 
Subject headings: stars: activity — stars: coronae — stars: individual (YY Men) — stars:flare — stars: late-type — 
X-rays: stars 



1. INTRODUCTION 

Magnetic activity is ubiquitous in late-type stars, although 
its level can vary dramatically. A common activity indica- 
tor is the ratio of the X-ray luminosity to the stellar bolomet- 
ric luminosity, Lx/i'boi; it varies typically from 10"^ in in- 
active stars to lO"-'- 10"^ in the most active stars (see, e.g., 
Ipav ata & Micela 2003 and Gudel 2004 for recent reviews on 
the X-ray emission of stellar coronae). Giant stars present 
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a peculiar behavior in which magnetic activity is strong for 
spectral types earlier than typically KO-2 III, while it van- 
ishes rapidly for cooler spectral types (e.g., Linskv & Haisch] 
1979; Avresetal. 1981; Linskv 1985; Hunschetal. 19961 
and references therein). Rotational breaking occurs in gi- 
ants around spectral type GO III cGrav ,,1^89). Although 
the ori gin of the divide remains debated. iHiin sch & Schrode^ 
Jl996t suggested that the evolutionary history of magnetically 
active stars of different masses explains the X-ray dividing 
line (al so Schroder. Hiinsch. & Schmitt 1998). Alternatively, 
iRosner et aL C1995 ) proposed that the dividing line could be 
explained by the change from a predominantly closed mag- 
netic configuration in stars blueward of the dividing line to a 
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predominantly open configuration in cooler giants. The latter 
con figuration induces cool w i nds an d the absence of hot coro- 
nae. ' Holzwarth & Schiissled (120011) argued that flux tubes in 
giants of spectral types G7 to KO remain inside the stellar con- 
vection zone in a stable equilibrium. 

FK Comae stars form a loosely-defined class of rapidly 
rotating single G and K giant stars, whose outstanding 
property is a projected equatorial velocity measured up to 
160 km s"', in contrast to the expected maximum of 6 km s"' 
for giants. One of the leading theories to explain the extreme 
properties of FK Com stars suggests that they were formed 
by coalescence of a contact binary when one of the compo- 
nents entered into the giant stage ( Bodd & Rucinski"1981aL 
IBopd & Stencel 1981b; Rucinski 1990). Simon & DraQ 
0989) noted alternatively that dredge-up of angular mo- 
mentum during the early red giant phase could explain 
their rapid rotat ion. Magnetic activity in FK Com stars 
is very strong ("Bopp & Ste ncell Il981bt iFekel et all IT986t 
lRuttenlfT98 7: Simon & Drake 119891) . In flie X-ray regime, 
their X-ray luminosities are in the range from Lx = 10^*^ to 
10^' erg s"', i.e., Lx/^boi 10"^- 10"^ (e.g., Maado et ajj 
1990; Weltv & Ramsev 1994; de Meideros & Mavor 1995; 
Huenemoerder 1996; Hiinsch. Schmitt. & Vog es^ .1998; 
Gondoin 1999; Gondoin, Erd, & Lumb 2002). 

We present in this paper new high-resolution X-ray spec- 
tra and high signal-to-noise light curves of the FK Com-type 
star YY Mensae (HD 32918) obtained with the Chandra X- 
ray Observatory and its High-Energy Transmission Grating 
Spectrometer (HETGS), and with the XMM-Newton Observa- 
tory. With these observations, we aimed to compare the X-ray 
emission of this bright FK Com-type star with that of other 
magnetically active stars. We show that YY Men is among 
the stars with the hottest coronae, with a dominant plasma 
temperature around 20-40 MK. Furthermore, we investigated 
the elemental composition of YY Men's corona to study abun- 
dance anomalies. 

The paper is structured as follows: Section |2l describes 
the main properties of YY Men and previous observations; 
section |3l gives the observation details, whereas the data re- 
ducti on i s described in section l4l (Sect. ITTI for Chandra and 
Sect. l4.2l for XMM-Newton). In section|5] we provide an anal- 
ysis of the light curves. Section |6| describe s the procedures 
for our spectra l ana lysis for Chandra (Sect. 16. U and XMM- 
Newton (Sect. 16. 2> . We approached the spectral inversion 
problem with three different methods to obtain the emission 
measure distribution (EMD) and abundances in YY Men's 
corona (Sect. l63l . Section0includes a discussion of our re- 
sults: abundances ("Sect. fTTl . EMD (Sect. f7~2\ . line broad- 
ening (Sect. 17. 3> . electron densities (Sect. 17. 4> . and optical 
depth effects (Sect. ^3}. Finally, we give a summary and our 
conclusions in section|8l 

2. THE GIANT STAR YY MENSAE 

YY Mensae (Kl III p) is a strong Ca II emitter 

Shdelman & MacConnelll 119731) at a distance of 291 pc 
'errvman et al. 1997), with a mass of M = 2.2 M© (see 
ondoin 1999), and a radius o f j; = 8.8x lO" cm = \2.1Rp, 
based on Teff = 4, 700 K from Randich. Gratton. & Pallavicinil 
(fT99l and L h ni =2 .7 x 10^^ erg s"' = 70Lq from 
ICutispoto et alJ (Il992h . It belongs to the FK Com- 
type class of giant stars (!Collier Cameron 1982). It 
displays a photometric period of P = 9.55 days, a pro- 
iected rotational v elocity o f vsin; = 45 - 50 km s~' 
dCoUier CamerorJ Il982i: iPiskunov. Tuominen. & Viihul 



fT990t iGlebocki & Stawikowskil l20()(1). and an in- 
cUnation angle of ; = 65 ° (IPiskunov et alJ 119901) . 
iGrewing. Bianchi. & Cassatellal ( Il986h presented the 
lUE ultraviolet spectrum of YY Men, showing very strong 
ultraviolet (UV) emission lines originating from the chro- 
mosphere and the transition zone. They obtained YY Men's 
EMD in the range logr(K) = 3.8-5.3, emphasizing its 
strength in comparison to other magnetically active stars. 
Line ratios indicated characteristic electron densities in the 
transition zone of n,, = 3-4 x 10'" cm"^ dGrewing et alJ 
.1986) . Strong flares we re observed in the radio and optical, 
lasting for several davs (||^ l99 gt ^lj|l987: Buntonetal. 198^ 
ICutispoto. Pagano. & Rodonoin992) . A sti'ong P Cyg profile 
in the Ha line during a flare indicate d strong outflowin g 
material with a velocity of 230 km s~' (iBunton e t al. 1989*). 
Finally, spot coverage in YY Men was generally found in an 
equatorial belt (Piskunov et al. 199Qi), in c ontrast to th e large 
polar caps found in many active stars (e.g.,'Vogt^l988'). 

In the X-ray regime, Bedford, Elliott, & Evles ( 198^ re- 
ported X-rays from YY Men detected with the EXOSAT CMA 
instrument. Giidel et al. ( 1996) gave results of observations 
made with ROSAT and ASCA. They found indications for a 
hot (up to 3 keV) dominant coronal plasma. Flaring activ- 
ity was reported, although hampered by many interruptions 
due to the low orbits of the satellites. Abundance depletion 
in most elements was suggested from fits to the medium- 
resolution ASCA spectra. The X-ray luminosities varied be- 
tween the observations and ranged from logLx = 32.3 to 
32.65 erg s"' (Gudel et al. 1996). 

Our Chandra and XMM-Newton observations are generally 
consistent with the above results. Nevertheless, deep, unin- 
terrupted monitoring and high spectral resolution provide the 
opportunity to significantly improve our understanding of the 
corona of YY Men. 

3. OBSERVATIONS 

Chandra and XMM-Newton observed YY Men as part of 
Cycle 4 and of the Reflection Grating Spectrometer (RGS) 
Guaranteed Time, respectively. We provide a log of the ob- 
servations in Tabled 

Chandra observed the giant with the Advanced CCD Imag- 
ing Specti'omete r (ACIS) in its 1x6 arr ay (ACIS-S) with the 
HETG inserted ( Weisskopf et al. 2002). This configuration 
provides high-resolution X-ray spectra from 1.2 to 16 A for 
the High-Energy Grating (HEG) spectrum and from 2.5 to 
3 1 A for the Medium-Energy Grating (MEG) spectrum. The 
HETGS has constant resolution (HEG: AA - 12 mA, MEG: 
AA ^ 23 mA FWHM). Thanks to the energy resolution of the 
ACIS camera, spectral orders can be separated. However, be- 
cause of the drop in effective area with higher spectral orders, 
we used only the first-order spectra. Further details on the 
instruments can be found in the Chandra Proposers' Obser- 
vatory Guide ' . 

XMM-Newton dJansen et al J 120011) observed YY Men in 
the X-rays with the RGS (A ~ 6-38 A with AA - 60- 
76 mA; den Herder et al. 2001) and the European Photon 
Imaging Camera (EPIC) MOS only (0.15-10 keV with 
£'/A£' = 20-50; Turner et al. 2001), since the EPIC pn cam- 
era ( St riider et al . 2001) was off-line due to a sudden switch- 
ing off of one quadrant (M. Guainazzi, priv. comm.). Simul- 
taneous optical coverage was obtained with the Optical Mon- 

* |http : //cxc .harvard. edu/proposer/POG I 
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itor (OM; iMason et al.ll200lh . although we did not make use 
of the data because YY Men was too bright (V = 8.1) for the 
U filter. Furthermore, we note that, due to a small pointing 
misahgnment, YY Men fell off the small 11" x 11 "5 timing 
window. 

4. DATA REDUCTION 
4. 1 . The Chandra data 

The Chandra data were reduced with the Chandra Interac- 
tive Analysis of Observations (CIAO) version 2.3 in conjunc- 
tion with the calibration database (CALDB) 2.21. We started 
from the level 1 event file. We used standard techniques de- 
scribed in analysis threads^. In particular, we applied a cor- 
rection for charge transfer inefficiency, applied pulse height 
analysis and pixel randomizations, and we destreaked CCD 
8 (ACIS-S4). Furthermore, we used non-default masks and 
background spectral extraction masks to reflect the modifica- 
tions introduced in the recent CIAO 3 release^. 

We then calculated grating response matrix files (RMFs) 
and grating ancillary response files (ARFs). CALDB 2.21 
contained separate line spread functions (LSF) for the pos- 
itive and the negative order spectra. We then co-added the 
positive and negative order spectra to increase statistics (using 
the co-added grating ARFs), and further binned the spectra by 
a factor of two. Our final HEG and MEG first-order spectra 
have bin widths of 2.5 and 5 mA, respectively. The effective 
exposure was 74.2 ks. Finally, we note that no contaminating 
X-ray source was detected in the zeroth order image. 

4.2. The XMM-Newton data 

The XMM-Newton X-ray data were reduced with the Sci- 
ence Analysis System (SAS) version 5.4.1 using calibration 
files from August 2003. The RGS data were reduced with 
rgsproc. The source spatial mask included 95% of the cross- 
dispersion function, whereas the background spatial mask was 
taken above and below the source, by excluding 97% of the 
cross-dispersion function. The mask in the dispersion-CCD 
energy space selected first-order events only and included 
95% of the pulse-invariant energy distribution. We generated 
RGS RMFs with 6,000 energy bins. The maximum effective 
exposure for the RGS was about 85 ks. 

Among the EPIC data, we used the EPIC MOSl data only 
since we preferred to give maximum weight to the high- 
resolution RGS spectra. We removed short periods of solar 
flare activity (4.25 ks), leaving ^ 81 ks of MOSl exposure 
time. We extracted the MOSl data from a circle (radius of 
47") around the source, and the background from a source- 
free region of the same size on a nearby outer CCD. We cor- 
rected the exposure for vignetting in the background spec- 
trum. We finally created MOSl RMF and ARE files for the 
source, using default parameters in SAS 5.4.1 (except for the 
detector map type, for which we used the observed instrumen- 
tal PSF at 3 keV). 

5. LIGHT CURVES 

We extracted light curves from the Chandra and the XMM- 
Newton observations. Since the zeroth order Chandra light 
curve was piled up, we used the dispersed first-order photons 
instead. The background for the dispersed photons was taken 



|http : //cxc .harvard. ed u/ciao/threads 

^ We used width_f actor_hetg=35 with task tg_create_mask, and lat- 
ter used max_upbkg_tg_d=6 .Oe-3 with task tgextract. 



from two rectangles "above" and "below" the source in the 
dispersion/cross-dispersion space, totaling an area eight times 
larger than the source. The scaled Chandra background was, 
however, about 150 times fainter than the source. The EPIC 
MOS 1 light curve is based on the extraction regions described 
in Sect. 14.21 Note that the light curves are not corrected for 
deadtime which is in any case negligible. 

Figure [2 shows the XMM-Newton MOSl, and the Chan- 
dra MEG and HEG first-order light curves with a time bin 
size of 500 s. YY Men was about 2 times brighter (logLx ~ 
32.5 erg s"') during the Chandra observation than during the 
XMM-Newton observation (logLx ~ 32.2 erg s"'), as derived 
from fits to the average spectra (Sect. |6jl. Both the XMM- 
Newton and Chandra light curves show no obvious flare, but 
they display a slow decrease in flux with time. The modula- 
tion is weak (15-20%), but a Kolmogorov-Smirnov test for 
the MOS data gives a very low probability of constant count 
rate (P <C 0.1). We also performed a Kolmogorov-Smirnov 
test to search for short-term variability in the XMM-Newton 
MOS light curve. None was found down to a time scale of 
- 300 s. 

The interpretation of the modulation is unclear. We believe 
that YY Men was not caught in the late phase of a flare decay, 
since there is no evidence of a temperature variation during 
the observations. However, it could be due to some rotational 
modulation effect in which active regions almost completely 
cover the surface of the very active YY Men. 

6. SPECTRAL ANALYSIS 

In this section we present the spectral analysis of YY Men's 
data taken with Chandra and XMM-Newton. Due to the domi- 
nant hot corona and to the high derived interstellar absorption 
(A^H ~ 7 X 10^" cm"2; Tab. 13 see below as well) toward YY 
Men, the Chandra HETGS spectrum is the best-suited grat- 
ing instrument; the XMM-Newton data were, however, most 
useful to access the long wavelength range (and thus the C 
abundance). 

6.1. Chandra 

We fitted the co-added first-order spectra of MEG and HEG 
simultaneously. However, we restricted the wavelength range 
to 1.6-16.2 A, and 1.8-25.0 A for HEG and MEG, respec- 
tively. Thanks to the low ACIS background during the obser- 
vation it was not necessary to subtract a background spectrum. 
Each background spectral bin contained typically 0-2 counts 
in HEG and — 4 in MEG /or an area 8 times larger than the 
source, thus the scaled contribution of the background is less 
than a half (a quarter) of a count per bin in MEG (HEG), much 
less than in the MEG or HEG first-order source spectra of YY 
Men. It allowed us to use the robust C statistics (Cash 1979) . 

6.2. XMM-Newton 

The XMM-Newton RGSl, RGS2, and MOSl data were fit- 
ted simultaneously"^. The spectra were grouped to contain a 

Although cross-calibration effects could, in principle, be of some impor- 
tance, we verified that it was not the case for YY Men by obtaining multi-T 
fits to the MOSl data alone, and to the RGS1+RGS2 spectra as well. The 
best-fit solution to the MOSl spectrum is close to the solution reported in 
Table |2| for the combined spectra. On the other hand, due to the high T 
of YY Men, the RGS data alone were not sufficient to constrain adequately 
the high-T component, and consequently absolute abundances are different. 
Nevertheless, abundance ratios relative to Fe are similar to the ratios found for 
the combined fit. We therefore feel confident that the combined RGS+MOS 
fits are not sfi'ongly biased by cross-calibration effects. 
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minimum of 25 counts per bin. We used the RGS spectra 
longward of 8.3 A, and discarded the EPIC data longward 
of 15 A due to the lack of spectral resolution of the CCD 
spectrum; in addition, we preferred to put more weight on 
the high-resolution RGS data. As the XMM-Newton data are 
background-subtracted, we used the statistics. 

6.3. Methodology 

We approached the spectral analysis with three different 
methods to obtain the EMD and abundances in YY Men's 
corona. We aim to compare the different outputs in order to 
discuss the robustness of the results. For each method, the 
coronal abundances in YY Men were compared to the sola r 
photospheric abundances given in lGrevesse & SauvalHl998l) . 
The first method uses the classical multi-T model to discretize 
the EMD. The second method obtains a continuous EMD de- 
scribed by Chebychev polynomials. The third method uses 
selected emission line fluxes to derive a continuous EMD and 
abundances. The continuum contribution is estimated from 
line-free continuum spectral bins. We describe the above 
methods in the following sections in detail. 

6.3.1. Method One: A Multi-Temperature Discretization of 
the Emission Measure Distribution 

This method uses the classical approach of a multi-T com- 
ponent model as a discretization of the EMD. The model, 
however, does not provide a physically reasonable description 
of the corona of a star since EMDs are thought to be contin- 
uous. Nevertheless, this approach is generally sufficient to 
obtain abundances (e.g., Schmitt & Ness 2004). 

We applied method 1 to the Chandra and XMM-Newton 
data. We used the Interactive Spectral Interpretation Sys- 
tem (ISIS) software versio n 1.1.3 ( Houck & Denicola 2000) 
for Chandra and XSPEC (lArnaudlll996 ) for XMM-Newton. 
The collisional ionization equilibrium models were based on 
the Astrophysical Plasma Emission Code (APEC) version 
1.3 jSmith et al. 2001), which was available in both soft- 
ware packages. Coronal abundances were left free, in addi- 
tion to the T and EM. We included a photoelectric absorp- 
tion component left free to vary that uses cross-sections from 
iMorrison & McCammonl (II983l) . We included thermal broad- 
ening for the Chandra spectra. While such a broadening is 
negligibl e for most lines, it plays some role for certain ions 
(see Sect. 1731. 

Figures l2aH2cl show the Chandra HEG and MEG spectra 
with the best-fit 3-r model overlaid. The overall agreement is 
very good. Figure|3]shows the XMM-Newton spectra with the 
best-fit A-T model overlaid. The higher sensitivity of XMM- 
Newton RGS at long wavelengths over Chandra HETGS al- 
lowed us to detect a low-T 3 MK) plasma component 
which was not detected with Chandra. Although the com- 
ponent is faint (as expected from the Chandra data), its inclu- 
sion decreases the statistics significantly (Ax^ = 32 for 2 
additional degrees of freedom). 

6.3.2. Method Two: A Continuous Emission Measure 
Distribution from Chebychev Polynomials 

This method obtains a continuous EMD described by 
Chebychev polynomials. An additional constraint for the 
spectral inversion problem is to keep EMs positive. This was 
achieved by approximating the EMD with the exponential of 
a polynomial as described by Lemen et al. ( .1989) . We use 
the convention that the differential emission measure, ip{T), 



is defined as 

<p(7') = «H«e^ (cm-^K-'). (1) 
dT 

Thus the total EM is given by EMjot = J ip(T)dT = 
J (p{T)Td{lnT). A graphical representation of the EMD is 
therefore given by 

EMD(r) = (^(r)rAlogr (cm-^). (2) 

In this paper, we preferred to plot EMD(r)/Alogr instead, 
since it is independent of the grid bin size (AlogT), to help 
comparison between the EMD reconstructed with methods 
2 and 3. We constructed a local model in XSPEC (Arnai^ 
[l996) in which the maximum polynomial degree can be fixed, 
and in which the T grid can be given as well^. Coronal abun- 
dances and photoelectric absorption were also left free to vary. 
Our model uses the same APEC version as in ISIS. We applied 
method 2 to both the Chandra and XMM-Newton spectra. 

We report here the results obtained with a grid between 
log r = 6.0 and 7.95 K and a bin width of A log T = 0. 15 dex, 
and for a polynomial degree of « = 7. Various combinations of 
T ranges, grids, and polynomial degrees were tested as well 
but they provided no improvement of the statistics. The T grid 
is fixed, and EMs are obtained from Chebychev coefficients a^ 
(k= 1 , . . . , 7) and from a normalization factor in XSPEC. We 
determined uncertainties on EMs for each T bin as follows: 
we obtained 68% (AC or Ax^ = 1) confidence ranges for the 
Chebychev polynomials ai , . . . , a?, and assumed that each co- 
efficient was normally distributed with a mean, /i,, equal to 
the best-fit value and with a standard deviation, cr,, equal to 
the geometric mean of its 68% uncertainties. We then used a 
Monte-Carlo approach and, for each coefficient, we generated 
10,000 pseudo-random values following a normal distribu- 
tion M{pi,cFi) (since Chebychev coefficients are in the range 
[-1,1], we assigned a value of -1 if the randomly generated 
value was < -1, and assigned +1 if the random number was 
> 1). Our Monte-Carlo approach provided us with 10,000 val- 
ues of the EM per T bin. We investigated the distribution of 
each EM per bin and noted that its logarithmic value closely 
followed a normal distribution. We therefore fitted the 13 dis- 
tributions with normal functions, and associated the derived 
standard deviation with the uncertainty of the EM per bin. 

6.3.3. Method Three: A Continuous Emission Measure 
Distribution Derived from Fe Emission Lines 

The last method used extracted fluxes from specific, bright 
lines to reconstruct the EMD. We applied it only to the Chan- 
dra data given the low line-to-continuum ratios and the broad 
wings in the XMM-Newton RGS spectrum, both of which 
make a clean flux extraction challenging in this case. 

To extract line fluxes, we used the Chandra HEGl and 
MEGl data simultaneously to fit each emission line to- 
gether with a "line-free" continuum in ISIS. The contin- 
uum model consisted of a single T , absorbed model con- 
taining the two-photon emission, radiative recombination, 
and bremsstrahlung continua and a pseudo-continuum. The 
pseudo-continuum in APEC consists of lines which are too 
weak to list individually and whose contributions are stored 
as a continuum. The fitting procedure consistently found a 

XSPEC has similar models, e.g., c6pvmkl, however, the models use the 
MEKAL database, a fixed T grid from log T = 5.5 to 8.0 with A logT = 0.10, 
and a sixth-order polynomial only. 
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best-fit logr = 7.53 K with log EM = 55.25 ±0.02 cm"^ and 
an absorption column density of logA^H = 20.89 ± 0.02 cm"^. 

Next, we applied an iterative EMD reconstruction method 
using the above line fluxes and the same APEC database 
as with the other methods. To allow for comparison with 
method 2, we used EMD(r)/Alogr for our graphical rep- 
re sentation of the EMD . The method is described in detail 
in lTelleschi et alJ ( 12004ft : however, we briefly summarize the 
principal steps. We treated unresolved line blends containing 
essentially only lines of one element like single lines, i.e., we 
computed new T-dependent emissivities for each considered 
line blend. Our EMD reconstruction starts by considering the 
Fe lines of Fe line blends. An approximate, smooth EMD 
estimated from the emissivities at the maximum line forma- 
tion T for each line serves as an initial approximation to the 
solution. The coolest portion of the EMD is estimated from 
the flux ratio of the O VIII to the O VII resonance lines. The 
EMD is binned into bins of width AlogT = 0.1 dex in the T 
range between log r(K) = 6-8. Once an EMD is defined, line 
fluxes were predicted by integrating the emissivities across the 
EMD. We then iteratively corrected the EM in each bin, using 
the iteration algorithm described by Withbroe ( 1975). The it- 
eration is terminated once the (for the deviations between 
predicted and measured fluxes) is no longer significantly im- 
proving, or if the reduced is £ 1 • At this point, the EMD 
has been determined only up to a normalization factor that 
directly depends on the absolute [Fe/H] abundance. To deter- 
mine these two quantities, we computed the predicted spectra 
for various values of [Fe/H] until the predicted continuum was 
in close agreement with the observations. We repeated this 
entire analysis twenty times, each time perturbing the mea- 
sured line fluxes according to their measurement errors and an 
assumed systematic uncertainty in the line emissivities (10% 
for each line). The scatter in the solutions provided an esti- 
mate of the uncertainties in the EMD. Abundance uncertain- 
ties include i) the measurement errors of the fluxes of the line 
blends, ii) the scatter in the abundance values derived from the 
different lines of the same element, iii) the scatter in the abun- 
dance values derived from the twenty EMD realizations, and, 
iv) for the absolute Fe abundance, the statistical uncertainty 
in the adjustment to the observed continuum level. 

6.3.4. Systematic Uncertainties 

Table |2l shows the best-fit parameters for Chandra and 
XMM-Newton for the 3 methods discussed above. We provide 
for each fit parameter an estimate of the uncertainty based on 
the confidence ranges for a single parameter of interest (90 %, 
except for method 3 where we quote l-cr ranges). However, 
we emphasize that such uncertainties are purely statistical, 
and do not inc lude systematic uncertainties except for method 
3, see ^6.3.31 Those are typically of instrumental nature (e.g., 
cross-calibration of the MEG and HEG, wavelength scale, ef- 
fective area). In addition, models using atomic data do not 
usually include uncertainties for the atomic parameters (e.g., 
transition wavelengths, coUisional rates, etc). These uncer- 
tainties (of the order of 10-50%; Lamina 2002) vary from ele- 
ment to element and from transition to transition, and thus are 
difficult to estimate as a whole. Certain groups of atomic tran- 
sitions, like low-Z L-shell transition s, are also simply m iss- 
ing in many atomic databases (e.g.. 'Lepso n et al]l2003h . If 
the temperature structure is such that those transitions are 
bright enough to be measured in an X-ray spectrum, their ab- 
sence in the atomic database can have a significant impact on 
the determination of the EMD and of abundances. Typically 



those lines are maximally formed at cool temperatures and 
can be found longward of 25 A (e.g., lAudard et al.ll200ral) . 
thus they are weak for most hot coronal sources (an example 
where weak lines are important is the cool F subgiant Pro- 
cyon; Raa ssen et al.l l2002). The Chandra and XMM-Newton 
X-ray spectra of YY Men's hot corona are, therefore, not 
strongly contaminated by these L-shell lines, which allowed 
us to include the relevant wavelength ranges. A small contam- 
ination could potentially still contribute to the systematic un- 
certainty of coronal abundances. In addition, since the spec- 
tral inversion problem to construct an EMD is mathematically 
ill-posed, EMDs are not unique, a nd several realizations ca n 
reproduce the observed spectrum (* Craig & Brown >'1976a'bl). 
Finally, uncertainties in the solar photospheric abundances ex- 
ist as well. All the above contribute to systematic uncertain- 
ties. Consequently a systematic uncertainty of at least 0.1 dex 
in the EMs and abundances should be expected. 

7. DISCUSSION 

In this section, we discuss quantities derived from the 
spectral and photometric analyses with Chandra and XMM- 
Newton. We provide comparisons between the observations 
obtained at two different epochs when relevant. 

7.1. Abundances 

Our analyses of the Chandra HETGS and XMM-Newton 
spectra showed a marked depletion of metals in YY Men's 
corona when compared to the solar photospheric composition 
JGrevesse & Sauvailll998h . There is generally good agree- 
ment between the abundances obtained with the various meth- 
ods. Below we compare the abundances obtained in this paper 
and discuss them in the context of an abundance pattern in YY 
Men's corona. 

7.1.1. Systematic Effects 

Various methods applied to the same data set gave relatively 
similar abundances within the error bars; however, small sys- 
tematic offsets can be observed. For example, abundances 
with Chandra data from method 2 are slightly lower with re- 
spect to those obtained from method 1 (Fig.|4t). This is prob- 
ably a trade-off of the fitting procedure between the EMD dis- 
cretization and the values of the abundances. On the other 
hand, no significant shift is visible in the XMM-Newton data 
(Fig.S)). Such a comparison is useful, since it shows that de- 
spite using the same data set, the energy range, and atomic 
database, systematic offs ets in ab undances can be derived 
from spectral fits (Craig & Brownii 1976b,) . We emphasize, 
however, that such an offset is small in the case of YY Men 
where we find excellent agreement. 

Method 3 was applied to the Chandra data only, and by 
definition, abundance ratios with respect to Fe are obtained, 
and later the absolute Fe abundance is determined from the 
continuum. We compare the abundance ratios obtained with 
methods 1 and 3 in Figure |^ and with methods 2 and 3 
in Figure Excellent agreement is derived for most ele- 
ments, demonstrating that abundance ratios are very stable 
regardless of the method used. We find no preference for ei- 
ther (so-called) global fits or line-based analysis. We note, 
nevertheless, weaker agreement for the Ca abundance ratio, 
which is larger with method 3 than with both methods 1 and 
2. The reported abundances with method 3 are weighted av- 
erages determined from the available lines of an element, as 
derived from the line fluxes and the shape of the EMD. Ca 
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was, however, determined from Ca XIX only because Ca XX 
is not obvious in the spectra. On the other hand, the first two 
methods could underestimate the Ca abundance. Indeed, the 
Ca XIX triplet is weaker in the model than in the MEG spec- 
trum (Fig.l2al. The robustness of abundance ratios, even with 
the simplistic method 1, is an important result that can be ex- 
plained by the fact that our models (including the multi-T ap- 
proach) sampled the EMD adequately to cover the cooling 
function of each line. 

We also compared the abundances obtained with the same 
method and atomic database, but with different data sets 
(Fig-13- The top panels in Figure|5]show the absolute abun- 
dances (relative to H) from the Chandra data on the x-axis and 
from the XMM-Newton data on the y-axis for the first method 
(left) and the second method (right). The XMM-Newton ab- 
solute abundances are generally lower by ^ 0.15 dex with re- 
spect to the Chandra abundances, regardless of the method. 
In contrast, we find a much better agreement when abundance 
ratios (e.g., with respect to Fe; bottom panels) are used. It 
again indicates that such ratios are significantly more stable. 

Although the absolute N abundance is larger with Chan- 
dra than with XMM-Newton, we find that the N/Fe abun- 
dance ratios match well within the uncertainties. Neverthe- 
less, with Chandra data, method 3 found a closer agreement 
to the XMM-Newton ratio than methods 1 and 2, for which the 
best-fit ratios are slightly larger. We attribute this effect to the 
faintness of the O VII triplet in the Chandra data, and the lack 
of coverage at long wavelengths. Both methods 1 and 2 could 
not constrain accurately the EMD at low temperature (where 
N VII is maximally formed), and thus the N abundance. This 
is, indeed, reflected in the lack of cool (^ 3 MK) plasma in the 
spectral fit to Chandra with method 1 . This component was 
required with the XMM-Newton data and thus the fit needed a 
lower N abundance to describe the N VII Lya line. Since the 
algorithm for method 3 includes the ratio of the O VIII line to 
the O VII line, there is a better coverage of the EMD at low 
temperature. 

The abundances of Ar and S are lower with XMM-Newton 
than with Chandra. Since lines are found mostly from K- 
shell transitions in the EPIC CCD spectrum, possible cross- 
calibration errors between the RGS and the EPIC could ex- 
plain the discrepancy. We note, however, that the upper limit 
of the Ca abundance with XMM-Newton is consistent with the 
abundance obtained with Chandra. We note that an alterna- 
tive solution could be that true abundance variations occurred 
in such elements between the XMM-Newton and Chandra ob- 
servations. Although we cannot strictly discard such an ex- 
planation, we find it improbable. 

7.1.2. An Inverse FIP Effect 

In contrast to the solar First Ionization Potential (FIP) 
effect in which low-FIP (< 10 eV) elements are overabun- 
dant with respect to the solar photospheric composition 
and in which h igh-FIP e l ement s are of photospheric com- 
position (e.g., Feldman '1992'; ' Laming. Drake. & Widing I 
[^95; Feldman & Lamins 2000), recent analyses of grat- 
ing spectroscopic data with XMM-Newton and Chandra 
have emphasized the presence of an inv erse d istribution in 
very a ctive stars (e.g., i Brinkman et al] 1200 It [D rake et al. 
200?, 'Audard, Gudel, & Mewe 2001b; "Gudel et al. 2001; 
Huenemoerder. Canizares. & Schulz 2001; Audard et al. 
2003i) in which the low-FIP elements are depleted with 
respect to the high-FIP elements. Unknown or uncer- 
tain photospheric abundances are, however, a general 



problem in magnetically active stars. Possibly, coronal 
abundances relative to the stellar phot ospheric abundances 
could show different patterns (e^., lAudard et al.l l2003t 
Sanz-Forcada, Favata, & Micela 2003). However, a study of 
solar analogs with solar photospheric composition, and of 
various ages and activity levels showed a transition from a 
solar-lik e FIP effect to an inverse FIP effect with increasing 
activity JGiidel et al.ll2002i) . Bright, active RS CVn binaries 
appear to follow the transition generally well ( Audard et al] 
l2003i) . Several studies have focused on patterns of abun- 
dances in magnetically active stars (e.g., Drake 1996). We 
refer the reader to recent reviews f or further information 
jDrakel2003atlFavata & Micelal2003blGiideL2004l. 

Figure plots the abundance ratios with respect to Fe in 
YY Men's corona, relative to the solar photospheric com- 
position (Grevesse & Sauval 1998). As shown previously, 
abundance ratios are more robust than absolute abundances 
in spectral fits. The XMM-Newton data also give access to 
the C abundance, whereas the Chandra data provide an es- 
timate of the Al abundance. Three features can be observed 
immediately in Fig. |6j i) low-FIP elements (< 10-15 eV) 
show solar photospheric ratios, ii) high-FIP elements (essen- 
tially Ar and Ne, and possibly O) are overabundant with re- 
spect to solar ratios, iii) N is hig hly overabundant; this fea- 
ture is addressed in Sect. 17.1.31 The first two features are 
reminiscent of the inverse FIP effect seen in m any stars with 
high activity levels (e.g., Au dard et al.l l2003h . such as YY 
Men (Lx/ Lhn] y 10" ^ using Lboi = 2.7 x 10^^ erg s"' from 
Cutispoto e t alJll992l) . The low-FIP elemental abundance ra- 
tios are all consistent with a solar ratio. In contrast, several 
active stars show a broad U-shape in their FIP pattern, with 
Al and Ca slightly overabundant with respect to Fe, whereas 
the abunda nces incre ase gradu ally with increasing FIP (e.g., 
Huenemoer der et al J2 003 : Ost en etall20 03). 

Many studies of stellar abundances focused 
mainly on the [Fe/H] photospheric abundance 

(e.g., Pallavicini. Randich. & Giampajjgj [T22i 

[ Randich. Gr atton. & Pallavicini 1993; Cayrel de Strobel et alJ 
1 19971 ICavre l de Strobel, Soubiran, & Ralite 2001]). Conse- 
quently, coronal FIP patterns can often only be compared 
with the solar photospheric composition, except in (rare) 
cases where photospheric a bundances of some el ements are 
relativ ely well-kno wn (e.g., 'Drake. Laming. & Widing 1^2^ 
| 1997t iGiidel etalJl2'0 02: Raassen et al. 2002; Audai' d et alJ 
"2003* ISanz-Forcada. Favata. & Micela .2004 .TeUeschi et alJ 
2004:)- YY Men is no exception. To our knowledge, 
only th e [Fe/H] abundance is available in the literature. 
iRandich et al. (1993.) quote [Fe/H] = -0.5; however, in their 
spectrum synthesis analysis, they used a solar abundance 
of A(Fe) = 7.63, (S. Randich 2003, priv. comm.), while we 
use in this paper the Grevesse & Sauval ( 1998) abundance, 
A(Fe) = 7.50. Therefore, we wifl use [Fe/H] = -0.37 in YY 
Men's photosphere. Together with the absolute abundances, 
[Fe/H] ~ -0.5 to -0.65, in YY Men's corona from Table |2| 
we conclude that YY Men's corona is depleted in Fe with 
respect to its stellar photosphere by about -0.15 to -0.30 dex. 
Assuming that the abundances of the other elements are 
depleted by a similar amount (-0.37 dex) in YY Men's 
photosphere, we cautiously suggest that low-FIP elements 
are depleted in its corona relative to its photosphere as well, 
whereas high-FIP elements are either photospheric or slightly 
enhanced. Obviously, additional photospheric data are 
needed to substantiate this claim; in particular, photospheric 
abundances from C, N, and O would be most useful. 
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7.1.3. The CNO Cycle 

Although most elements seem to follow an inverse FIP ef- 
fect pattern, another explanation is required for the high abun- 
dance of N (absolute or relative to Fe) in YY Men. Indeed, the 
N/Fe abundance ratio in Fig.|6lis much larger than the O/Fe 
abundance ratio, despite similar FIPs for O and N. Some en- 
hancement could be due to the inverse FIP effect, but we es- 
timate this effect to no more than ^ 0.15 dex, based on the 
[O/Fe] ratio. Nitrogen enrichment in the photosphere of gi- 
ant stars due to dredge-up in the giant phase of material pro- 
cesse d by t he CN O cyc le in the stellar interior is a candidate 
(e.g., EES [HM 

Il967ft . Significant mass transfer in close 
binaries can also reveal the composition of the stellar inte- 
rior Such stars show X-ray spectra with an enhanced N/C 
abundance ratio due to the CNO cycle ( Schmitt & Ness 2002; 
iDrake 2003a b; Schmitt & Ness 2004). Similar enhancement 
in V471 Tau was reported by Drake & Sarna (>2003c) as ev- 
idence of accreted material during the common enveloppe 
phase of the binary system. A high N/C ratio has been found 
in hot stars as well (e.g., Kahn et al. 2001; Mewe et al. 2003}). 
YY Men, however, is a single giant, although it could have 
been formed in the merging of a contact binary when one 
component moved into the giant phase (e.g., Rucinski 1990) 

The abundance ratio is [N/C] +0.70;!:" ^ based on the 
XMM-Newton fits reported in TableQ Schaller et al. ( 1992) 
report an enhancement of [N/C] =+0.56 for a star of ini- 
tial mass M = 2.5 Mq, and metalHcity Z = 0.02 {X = 0.68, 
Y = 0.30). Due to the uncertain origin of YY Men, such 
a number is indicative only; however, it suggests that the 
high N coronal abundance in fact reflects the photospheric 
composition of YY Men. Interestingly, the prototype of its 
class, FK Com (G5 111), shows an unmixed N V/C IV ratio 
(Ifekel & Balachandran 19931). This is consisten t with the line 
ratio of N VII/C VI in X-rays (Gondoin et al.''2002). These 
ratios point to a possible evolutionary difference between 
FK Com and YY Men, in which the former is still in the 
Hertzsprung gap phase and has not undergone the He flash, 
and YY Men is a post-He flash clump giant with strong X-ray 
emission . 

1 ISchaller et all Jl992t . the flieoretical C/Ne ra- 
tio after the first dredge-up of a giant is also depleted 
(-0.16 dex), whereas the O/Ne ratio essentially remains 
constant (-0.02 dex). From our fits, [C/Ne] ~ -0.5, and 
[O/Ne] ~ -0.5, lower than expected from evolutionary codes, 
but with large error bars. Nevertheless, we recall the enrich- 
ment of coronal Ne by a factor of about +0.3 to +0.5 de x in 
active stars f^g.- lBrinkman et alJ200U|Drake et alJ2001l ) due 
to the inverse FIP effect. Such an enhancement lowers signif- 
icantly the abundance ratios relative to Ne, which can explain 
the lower abundance ratios of [C/Ne] and [O/Ne] than the evo- 
lutionary predictions. 

The CNO cycle keeps the total number of C, N, and 
O atoms constant. In addition, for a star of intermediate 
mass (2.5 Mp,), the number of Ne atoms does not change 
(^challer et al. 1992). Thus, for a star with an initial so- 
lar composition, the ratio (Ac + An + Ao)/ANe should re- 
main solar, where Ax is the abundance number ratio of el- 
ement X relative to H. The solar photospheric values of 
iGrevesse & SauvalUi^gsl) are Ar = 3.31 x 10"^ An = 8.32 x 
10"^ Ao = 6.76 X 10"^ and Ane = 1 .20 x 10"^ yielding a ra- 
tio of 9.1. Our best-fit abundances from XMM-Newton data 
yield a ratio of ~ 4.3 with an estimated uncertainty of ±2.5. 
Although the lower ratio again could cast some doubt on the 



detection of CNO cycle in YY Men, we reiterate that coronal 
Ne is enhanced due to the inverse FIP effect. If the photo- 
spheric Ne abundance in YY Men is two times lower than in 
its corona, the above ratio would be in good agreement with 
the solar ratio. Finally, although the above discussion assumes 
near-solar abundances in the stellar photosphere , the o bserved 
Fe depletion in the photosphere ( Randich et aO l 19931) argues 
against a solar composition in YY Men in general. 

7.2. Emission Measure Distribution 

Figure displays the EMD of YY Men's corona derived 
from our spectral fits to the Chandra and XMM-Newton data. 
We reiterate that we use EMD(r)/Alogr for our graphical 
representation of the EMD for methods 2 and 3, whereas we 
plot the EM of each isothermal component for method 1 . Al- 
though the plotted EMs for the latter method lie below the 
EMs per bin curves for the two other methods, we empha- 
size that this is solely for plotting purposes since the total 
EMs in YY Men are similar with the three method (logEMtOt 
(cm"^) 55.3 and 55.1 during the Chandra and XMM- 
Newton observations, respectively). The upper panel with a 
linear ordinate emphasizes the dominant very hot plasma in 
YY Men. The lower panel, using a logarithmic scale, reveals 
the weak EM at lower temperatures (< 20 MK). Other FK 
Com-type stars show similar EMDs with a dominant plasma 
at r > 20 MK (Gondoin et al. 2002; Gondoin 2004). 

The various methods show consistently an EMD peaking 
around 20 - 40 MK. However, for the same data set, the de- 
rived EMD can be slightly different. For example, the realiza- 
tion obtained from method 3 with the Chandra spectra show 
a smoother solution than the Chebychev solution (method 2). 
Method 3 starts with a smooth approximation and introduces 
structure only as far as need to fit the line fluxes. A suffi- 
ciently good was reached before any strong modification 
of the EMD at high-T was introduced. The good statistics 
and the excellent match between the abundance ratios show 
that both EMDs are realistic and valid. We see no preference 
for either method, in particular in the light of the considerable 
systematic uncertainties and the consequent ill-posed spectral 
inversion problem. 

We also note that the dominant plasma is at slightly lower 
temperatures (20 MK) with XMM-Newton than with Chan- 
dra (40 MK). A larger EM for the highest-T component is 
possible with XMM-Newton but is unstable to the fitting al- 
gorithm. It remains unclear whether YY Men was actually 
cooler during the XMM-Newton observation. Indeed, the sta- 
tistical EM uncertainties with method 1 are smaller at high-T 
than at mid-T with Chandra but the trend is reversed with 
XMM-Newton (Tab. |2j, which suggests that, due to the spec- 
tral inversion problem, fits could have preferred a somewhat 
cooler realization for the EMD of YY Men during the XMM- 
Newton observation. 

7.2. 1 . Coronal Structures 

The lack of spatial resolution (in contrast to the Sun) 
does not allow us to understand what structures are present 
in the corona of YY Men. We are re stricted to lo op for- 
malisms (e.g., iRosner. Tucker. & Vaiana 1978; Mew^^lJ 
1982; Schriiver. Lemen. & Mewe 1989; van den Oor d et alJ 
1997). The analytical formulae allow us to derive some in- 
formation on the loop structures after placing certain assump- 
tions (e.g. constant pressure, similar maximum temperature, 
constant loop cross-section, etc). Discrepancies between the 
observed EMD and the one predicted by loop models allow 
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us to test some of the above assumptions and modify them 
accordingly if necessary (e.g., by allowing loop expansion 
towards the top or the presence of distinct families of loops 
with different maximum temperatures). The high coronal 
temperature (a few tens of MK) and the low surface gravity 
(~ 0.01^0) for YY Men imply that the pressure scale height 
Hp is large enough (^^ 5 x 10'^ cm = 67?*) to ensure constant 
pressure in most X-ray emitting loops. 

Argiroffi, Magaio, & Peres (2003|) and lScelsi et"an J2004 
have recently published EMDs for the RS CVn binary Capella 
observed by Chandra and for the single rapidly rotating 
GO III giant 31 Com observed by XMM-Newton, using the 
Markov-Chain Monte Carlo EMD reconstruction method by 
iKashvap & Drake ( 1998). Figure El compares the EMD de- 
rived by us for YY Men with the EMDs of Capella and 
31 Com derived by these authors. In the figure we have 
used the EMD derived from Chandra data using method 3 
(and the graphical representation given in Eq. |2}, since this 
method uses the sam e bi n width AlogT = 0.1 as used by 
lArgiroffi et aT ("2003") and lScelsi et aT (2004). Although the 
EMDs of Capella and 3 1 Com are similar, with 3 1 Com being 
only slightly hotter than Capella, the EMD of YY Men peaks 
at higher T and shows a shallow increase to the peak (Fig.^}. 

Figure |9] compares the EMD of YY Men derived from 
XMM-Newton data using method 1 with those derived with 
basically the same method (RGS + M0S2 data and APEC 
code) by Audard et al. (2003) for a sample of RS CVn bina- 
ries of various activity levels (including Capella). We have 
added for comparison t he res ults of a 3-7" fit to MOS2 data 
derived by IScelsi etakl (|2004') for 3 1 Com as well as the re- 
sults of 2-T and 3-T fits to EPIC data (MOS-npn) of two other 
FK Comae-typ e stars (V1 794 Cyg a nd FK Com itself) from 
ICiondoin et all ll2002h and iGondoinl ll2004 . Figures El and 
|9l emphasize the exceptional nature of YY Men (and of FK 
Com-type giants in general) compared to other active stars in 
terms of high T and high EM. 

We can use the approach proposed by Peres et al. ( 2001 ) to 
infer the properties of coronal loops from the observed EMD 
shape. If we assume the loop model of Rosner et al. ( 1978), 
we find that the shape of the emission measure distribution for 
a single loop does not depend on the length of the loop, but 
only on the loop maximum temperature T^Ti ax, and it s func- 
tional form is EMD(r) cx up to r^ax. P eres et alJ ll2001l) 
showed that, by grouping coronal loops by their maximum 
temperature irrespective of length, it is possible to describe 
the integrated properties of the corona, and the total emis- 
sion measure distribution, as the sum of the contributions of 
different families of loops, each characterized by a different 
value of the maximum temperature Tmax- They showed, in 
particular, that under these assumptions the ascending slope 
of the EMD(r) of the whole corona is linked to /3, whereas 
the power-law index of the descending slope is linked to the 
distribution of the maximum temperatures of different classes 
of coronal loops. 

For YY Men we find a flatter EMD increase (/3 ^1.5 from 
method 3 and ^ 2.15 for method 2; see Fig.0, closer to the 
solar slope. This could indicate that the heating is more uni- 
form in the loops of YY Men than in Capella and 3 1 Com 
(which have - 5-6; Argiroffi et al. 2003; Scelsi et al. 2004) 
and that the loop cross-section is approximately constant. We 
caution, however, that the method u sed to derive t he EMD for 
YY Men differ from those used by Argiroffi et al. (2003) and 
IScelsi et al . (2004). This could be responsible in part for the 
different slopes. A very steep gradient {T~^) is obtained at 



temperatures > Tmax for Capella and YY Men and somewhat 
le ss (but with l arger e rrors) for 3 1 Com. In the loop formalism 
of lPeres et alJ (|200'l*), such steep slopes imply that most of the 
loops have Tmax equal to, or only slightly larger than the peak 
value of the total emission measure distribution. 

7.2.2. Flare Heating 

The solar corona cools radiatively mostly around 1-2 MK, 
with very weak EM at higher T (e.g. .'Or lando. Peres. & Reala 
2000). Stellar coronae display higher coronal T, with 
FK Com-type stars at the upper end of the active stars. 
A major problem of the solar-stellar connection studies is 
to understand the mechanism at the origin of the high 
coronal T. Coronal heating by flares is an attractive 
mechanism because flares often display large tempe ratures 
(e.g., Giidel et al. 1999; Franciosini, Pallavicini, & Tag haferril 
i200L) . The concept assumes a statistical contribution of 
flares to the energy budget in solar and stellar coronae (e.g., 
'Parker"1988). Recent work on flaring stars has shown that 
a continual superposition of flares (which follow a power- 
law distribution in energy) could radiate sufficient energy 
to explain the observed X-ray luminosity, suggesting that 
flares are significant contrib utors to the coronal heating in 
the Sun and in stars (e.g., Giidel 1997; 'Krucker & Bend 
1998; Audard, Giidel, & Guinan 1999; Audard et al. 200(J 
Parnell & Jupp 2000; Kas hvap et al. 2002; Gudel et al. 2003). 

The radiative cooling curve, A(r), scaled to coronal abun- 
dances in YY Men, is flat in the 3-20 MK range ( A 1 .2 x 
10"^"' erg s~' cm-'), in contrast to the curve for solar abun- 
dances (Fig.llO>. The radiative loss time, = likT /{n^KiJ)), 
is thus essentially directly proportional to the temperature (as- 
suming constant «e, see Sect. I7.4> in this range. Figure 
shows the radiative loss time as a function of T, using the 
cooling curve A(r) in YY Men's corona, and assuming elec- 
tron densities of n^ = lO** cm"-' and «e = 10'" cm"-'. There is no 
evidence for h igh densities above 10^- lO"* cm"^ in YY Men 
(see Sect. l7.4t . The radiative loss time is long (> 0.1 days), 
thus we can expect that the occurrence time scale of stochastic 
flares is much shorter than the cooling time. 

Consequently, the flare heating mechanism appears attrac- 
tive in YY Men, with frequent hot flares replenishing the 
corona on time scales shorter than the radiative cooling time 
scale, thus keeping YY Men's corona hot by continual reheat- 
ing. The lack of obvious flares in the Chandra and XMM- 
Newton light curves over observing time spans of about a day, 
however, does not immediately support this hypothesis. The 
XMM-Newton EPIC MOSl light curve (Fig.[0 is on average 
/i ^ 1.3±0.1 ct s"', corresponding to Lx ^ 1.6 x 10^^ erg s"'. 
For a 3 CT detection of a flare, a minimum of 0.3 ct s"' is re- 
quired at peak, correspond to L^'^'^ > 4 x 10^' erg s"' at YY 
Men's distance. A similar reasoning applies to the Chandra 
light curve. Such flares are at the high end of the flare lumi- 
nosity scale in stellar coronae and about 1,000 times stronger 
than the strongest solar flares. Due to the large distance of YY 
Men, it is, therefore, no surprise that current detectors are un- 
able to detect small-to-moderate flares in its corona since they 
are overwhelmed by fluctuations due to Poisson statistics. 

The above result suggests that if the flare heating mecha- 
nism is operating in YY Men, the flare occurrence rate dis- 
tribution in energy, dN/dE oc £"", must be steep. GiidH 
( 1997) and Giidel et al. (2003) argued that the shape of the 
EMD can provide information on the power-law index a. In 
brief. iGudel et alJ ( l2003h obtained fliat the EMD is described 
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by two power laws, below and above a turnover. 



EMD(r) OC I j,_(,,_V)(o-2/3)/(l-/3)+2i-0^ 



T < Tq, 
T > To. 



(3) 



The turnover temperature, Tq, is related to the minimum 
peak luminosity, Ln,in = afT^~'^, from which the flare occur- 
rence rate distribution needs to be integrated. The EMD rises 
with an index 2/C. The index C corresponds to the index of the 
flare temperature decay as a function of the plasma density, 
T OC n '^. Typical values are C = 0.5-2 (Reale, Serio, & Peres 
Il993h . The other indices come from K(T) = fT-'^ (Fig. [Tot . 
the flare decay time t oc (E being the total radiative X- 
ray energy), and the flare peak emission measure, EM = aT'\ 
lAschwandenl (| 1999') reporte d fc = 7 and a = 0.6875 cm"^ K"^ 
for the Sun JFeldman et al] ll995t iFeldman 1996 quoted an 
exponential function, EM(r) = 1.7 x io" i37i+46 cm"^ which 
can be approximated with h ^ 5-1 in the range of interest; 
furthermore, Gudel 2004 obtai ns b = 4.3 zt 0.35 u sing stellar 
flares). We refer the reader to iGudel et alJ (l2003i^ for more 
details. 

In YY Men, Tq occurs around 30 MK, which corresponds 



to L„ 
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10^'' erg s-i (/ = 2.28 x lO'^^ erg s"' 



and = -1/2, = 7). Below Tq, EMD(r) cx implying 
C - 0.9. Above 7b, EMD(r) cx 7"^ For /3 = (flare decay 
time T independent of the flare energy), we obtain^ a ~ 2.7, 
weakly dependent of the exact choice of b (lower b values in- 
crease a). For a weak energy dependence, [3 = 0.25, we get 
a = 2.6, which still shows that a wealth of flares with small 
energies could dominate the coronal heating in YY Men, but 
remain undetected with the sensitivity of current detectors. It 
is worthwhile to note that L^i„ is of the order of very large so- 
lar flares, implying that, while smaller flares could be present 
in YY Men, they are not needed to explain the X-ray emission 
ofYY Men. 

7.3. Line Broadening 

Excess line width is observed with Chandra in Ne X 
A 1 2 . 1 3 A, and, albeit at lower significance, in Si XIV A6 . 1 8 A, 
O VIII A18.97 A, Mg XII A8.42 A, and possibly N VII 
A24.78 A. We give in Table |3] the line fluxes and the width 
measurements in excess of the instrumental width together 
with their 68% confidence ranges. Other emission lines, how- 
ever, do not show evidence of widths in excess of the instru- 
mental profile. We propose that the observed line broadening 
has a Doppler thermal broadening origin. In a Maxwellian ve- 
locity distribution of particles with a temperature T, thermal 
movements broaden the natural frequency (or wavelength) 
width of a transition c entered at vo, and produce a fine pro- 
file function (piv) (e.g.. lRvbicki & Lightmanll979h . 



-{v-vafl(AvDf 



where the Doppler thermal width is obtained from 



AAo 



2kT 



(4) 



(5) 



vq Ao 

with M as the atomic mass of the element. The observed 
broadening, a, is related to the Doppler thermal width, 

_ AAz5 _ Ao fi^ 
^~ ^ ~ c\ M' 



(6) 



* Here, we use the continuous EMD with method 2. Slightly lower indices 
a are found from the decreasing slope of the EMD with method 3; however, 
they still are above the critical value of 2. 



For an isothermal plasma of 40 MK, the Doppler broad- 
ening is 18.7 mA for C VI, 12.7 mA for N VII, 9.1 mA for 
O VIII, 5.2 mA for Ne X, 3.2 mA for Mg XII, 2.2 mA for 
Si XIV, 1.6 mA for S XVI, 1.1 mA for Ar XVIII, and 0.9 mA 
for Ca XX Lya lines. For Fe lines, the broadening is 2.6 mA 
and 3.8 mA at 10 A and 15 A, respectively. 

Although the above isothermal approximation works for the 
strongly peaked EMD of YY Men, a more thorough analysis 
can be given. Based on line emissivities and the derived EMD, 
we determine from which temperature the photon luminosity 
arises. The photon luminosity of line /, is determined as 



£, = / HT)^(T)dT 



(7) 



where A,(r) is the photon emissivity (ph cm^ s~') of line / 
from APEC, and tf{T) is as in Eq. [2 Therefore, a graphical 
representation of the photon luminosity distribution is Lff) = 
Ai(T)ip(T)T AlogT = Ai(T) x EMD(r), which is shown in 
Figure El (top panel) for O VIII and Ne X Lya lines, using 
abundances and EMD from the Chandra method 2 fit. De- 
spite the line maximum formation temperatures being at 3 MK 
and6 MK, respectively, most photons come from T bins in the 
range of 6 - 15 MK and 20-60 MK! Also, line photons in bin 
T are broadened by a specific square width, crfiT) (Eq. |6j. 
The distribution of square widths weighted by the photon lu- 
minosity, < (jf{T) >= (jfiT) X Li(T)/Ci, shows the contribu- 
tion of each T bin to the square width of the line (Fig. El 
bottom panel). The distribution shows that the width is domi- 
nated by the 30-50 MK plasma component. 

Since we used the Chandra calibration with the latest LSF 
profile, we believe that the measured excess widths are not 
related to inaccuracies in the instrumental profile. We tested 
further the robustness of our results, e.g., in Ne X and Si XIV. 
We simultaneously fitted the HEG negative and positive first- 
order line profiles with Gaussian functions, with the wave- 
lengths left free to vary in both spectra (to allow for wave- 
length calibration inaccuracies between the spectra), and let- 
ting the width and flux free as well, however linking the two 
parameters in both sides. Then, an instrumental profile was 
fitted to the emission line. The latter wavelengths and flux 
were used to simulate 10,000 realizations of a similar emis- 
sion line. Each realization was then fitted with a Gaussian 
function with free width. This procedure aims to obtain the 
fraction of best-fit Gaussian widths larger or equal to the de- 
tected excess width, which provides us with an estimate of the 
detection significance. 

The Ne X Lya line showed a width of a = 6.2 mA with a 
90% confidence range of [4.4,7.8] mA. The moments of the 
widths of the 10,000 Monte-Carlo simulations show an aver- 
age of 1.2 mA, with a standard deviation of 1.0 mA. More 
interestingly, the maximum width fitted from a model with an 
instrumental profile is 5.3 mA, lower than the best-fit excess 
width. Essentially, this shows that there is a probability less 
than 10"'* that the detected excess width in the Ne X Lya line 
is spurious. Taking into account the 90% confidence range, 
there are only 66 occurrences out of 10,000 where a best- 
fit width larger than or equal to the lower limit of the range, 
i.e., 4.4 mA was spuriously measured, although we empha- 
size that the corresponding line fluxes are weak. Figure [Ol 
(left) shows the AC = 1 and AC = 2.71 confidence maps in the 
line flux/sigma space together with the best-fit values of the 
10,000 Monte-Carlo simulations. Figure [Ol (right) shows the 
observed line profile, the instrumental profile, and the best-fit 
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Gaussian model. 

Similarly, the Si XIV Lya line showed some evidence of an 
excess width. However, the poor signal-to-noise ratio and the 
smaller spectral power at 6 A (essentially half the power at 
12 A) diminishes the significance of the detection. We ap- 
plied the same procedure as for the Ne X line; however, we 
rebinned the spectra by a factor of two (bin size of 5 mA) to 
increase the number of counts per bin. The best-fit width was 
4.4 mA with a 90% confidence range of [2.4,6.2] mA. The 
significance is lower (209 out of 10,000 values are larger than 
or equal to the best-fit width, and 2,107 out of 10,000 for the 
lower limit); however, the loci of the Gaussian widths of the 
Monte-Carlo realizations suggests that the detection may be 
real (Fig.fH] left panel). Figure [H] (right) again shows the 
observed line profile, the instrumental profile, and the best-fit 
Gaussian model. 

It appears that a broadening smaller than ^ 3 mA proves 
difficult to measure, which could explain the absence of 
broadening in the short wavelengths, e.g., in the S XV, 
Ar XVIII, Ca XX lines. On the other hand, it cannot explain 
the non-detections in, e.g., the Fe lines around 15 A and the 
upper limit of N VII at 24.7 A. We believe that the poor signal- 
to-noise ratio of these lines against the underlying continuum 
is at the origin of this discrepancy. Few counts are measured, 
in contrast to the strong signal of the Ne X Lya line. We note 
that, despite the strong flux of the O VIII Lya line, interstel- 
lar absorption and effective areas both conspired to reduce the 
amount of detected counts significantly. 

Although the above interpretation focused on Doppler ther- 
mal broadening only, we note that line broadening could be 
due to the stellar rotation as well, in part or completely. In- 
deed, Avres et al. ( 1998) found evidence in UV of broad line 
profiles of the fastest rotating gap giants (YY Men was not 
part of the sample, but it included FK Com). The broaden- 
ing suggested the presence of emission sources in the transi- 
tion zone at heights of ~ 7?* above the photosphere. In ad- 
dition, Chung et al. (2004) found excess line broadening in 
the Chandra X-ray spectrum of Algol which they interpreted 
as rotational broadening from a radially extended coronae at 
temperatures below 10 MK and with a scale height of order 
the stellar radius. 

Measurements of the projected equatorial rotational ve- 
locity in YY Men indicates that vsin/ = 45 km s"' 
dPiskunov et alJl99(l) . which implies a maximum wavelength 
shift of 

AA„ax=^vsin/=1.82x mA (8) 

with A in A. Therefore, for structures at the stellar surface, ro- 
tational broadening is smaller than Doppler broadening in the 
Chandra HETGS and XMM-Newton RGS wavelength range, 
and therefore it probably does not explain the observed line 
broadening. However, if coronal X-ray emitting material is 
high above the surface (e.g., at the pressure scale height. 
Hp = 6Ri,), it could produce significant broadening visible in 
the X-rays as well. If interpreted as purely rotational broad- 
ening, the excess width for Ne X, cr ^ 6 mA (Tab.O, implies 
a velo city of about 200 krn s~', using the inclination angle of 
65° dPiskunov et alj[T990l) . The latter authors also found no 
polar spots in YY Men but equatorial belts. Consequently, 
the excess width suggests that the coronal X-ray emitting ma- 
terial lies at the equator at heights of about 3/?* above the 
surface, i.e., below the pressure scale height. We note that at 



the wavelength of Si XIV, a larger velocity (300 km s"') is 
required to match the observed excess width, i.e., a higher al- 
titude (6/?*). Although rotational broadening remains a possi- 
ble explanation, we favor Doppler thermal broadening instead 
because of the dominant very hot plasma temperature in YY 
Men (see Fig. El bottom panel), in comparison, e.g., to the 
lower plasma temperature in Alg ol where rotational broaden- 
ing was the preferred scenario o f lChungetalJll2()0l . 

7.4. Densities 

The Chandra and XMM-Newton grating spectra cover 
ranges that include transitions whose intensities are density- 
sensitive. In particular, line ratios of the forbidden (/) to in- 
tercombination (/) lines of He-like transitions are most use- 
ful since they are sensitive to plasma electron densities cov- 
ering those found in stellar coronae (10^-10" cm"^; e.g., 
^ess et al. 2002, 2003a, 2004; Testa, Drake, & Peres 2004a). 

The very hot corona of YY Men is, however, problematic 
to derive electron densities from He-like transitions since the 
latter are faint and of low contrast against the strong underly- 
ing continuum (e.g., Fig. l2aH2ct . Indeed, most He-like triplets 
where the Chandra HETGS and XMM-Newton RGS effective 
areas and spectral powers are large enough (Si XIII, Mg XI, 
Ne IX, O VII) ai-e formed at T < 10 MK, where the EM is 
much lower (^ 1 dex) than the peak EM (Fig.Ql. In addition, 
blending is frequent (e.g., contamination of Ne IX by Fe lines; 
Ness et al. 2003b). 

We extracted individual line fluxes for Si XIII, Mg XI, 
Ne IX, and O VII from the Chandra spectrum (Tab. 13 because 
HETGS offers the best spectral resolution. We used the Chan- 
dra HEGl and MEGl data simultaneously to fit emission 
lines individually together with a "line-free" continuum (see 
Sect. 16.3.31 . No line fluxes of He-like triplets were obtained 
from the XMM-Newton RGS spectrum because of the lower 
spectral resolution for Si XIII and Mg XI, of the Fe blending 
in Ne IX, and of the low signal-to-noise ratio of the O VII and 
N VI triplets (Fig|3. 

Unfortunately, the lack of signal of the He-like triplets in 
YY Men produced large confidence ranges^. The derived 
R ratios (R = f/i) have uncertainties that do not allow us to 
constrain significantly the electron densities. However, the R 
ratios based on the best-fit fluxes give no indication of high 
(> 10'" cm"^) densities in YY Men. 

7.5. Opacities 

Our models assumed an optically thin plasma in YY Men's 
corona. A systematic study of stellar coronae with XMM- 
Newton and Chandra by Ness et al. (2003c) showed that stel- 
lar coronae are generally optically thin over a wide range of 
activity levels and average coronal temperatures. However, a 
recent study of the Lya/Ly /j line ratio in active RS CVn bina- 
ries bv lTesta et alJ ( I20041A found evidence of opacity effects, 
suggesting that opacity measurements by means of line ratios 
of Fe XVII could be hampered in active stellar coronae which 
are strongly Fe-depleted (e.g., , Audar d et aL 2003). 

We obtained estimates of the optical depths, t, from 
line ratios, using the "escape-factor" model designed by 

' We derived confidence ranges from a grid of line fluxes, allowing the 
statistics to reach AC = 1 . Our method accounts for uncertainties in the 
determination of the continuum, for Nh, and for calibration. We avoided 
calculating uncertainties from the square-root of the number of counts in the 
emission line since this method, though producing smaller uncertainties, does 
not take into account the effects mentioned above. 
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iKaastra & Mew3 (fT99l for r < 50, 



Rq l+0.43r' ' ' 

where R is the ratio of the flux of an opacity-sensitive line to 
the flux of an opacity-insensitive line, e.g., the ratio of Fe XVII 
A15.01 A to Fe XVII A15.26 A (Tab. and Ro is the ra- 
tio for an optically thin plasma. We note that the optically 
thin ratio of the above lines is subject to debate. Whereas 
theore tical codes range from 3.0-4.7 JBhatia & Doschel3 
1 1992V laboratory measurements obtain lower ratios (2.8-3.2; 
|Brown et al. 1998; Lamina et al. 2000). Despite these uncer- 
tainties, F(15.01)/F(15.26) ^ 2.7 is compatible with no sig- 
nificant optical depth in YY Men's corona (although uncer- 
tainties formally suggest a possible optical depth, however 
with r < 3). The ratio of the Fe XVII A 15.01 A to Fe XVII 
A16.78 A, a measure of optical depth as well, is close to 
the expecte d values f or an optically thin plasma dSmith et alJ 
ITOOl: Doron & BehaL.2002; Ness et eLL2003d) 

In view of the results bv iTesta et alJ J20()4bl) . opacity ef- 
fects could possibly be better measured in the Lya/Ly/3 ratios 
of abundant elements. As derived from Table |3 the line ratio 
for Ne X is consistent with the theoretical ratio in APEC. The 
ratio for O VIII is slightly lower, but we argue that contamina- 
tion by Fe XVIII could reduce the ratio artificially. Lya/Ly/3 
Une ratios of, e.g., N VII, Mg XII, Si XIV look consistent with 
the theoretical ratios as well, although the weakness of the 
Ly/3 lines makes the accurate measurement of such ratios dif- 
ficult. In conclusion, there is no strong support for a signifi- 
cant opt ical depth in the co rona of YY Men, in line with the 
study bv lNess et aLl(l2003cl) . 

8. SUMMARY AND CONCLUSIONS 

In this paper, we have presented our analysis of the X-ray 
emission of the FK Com-type giant star, YY Men, observed 
by Chandra HETGS and XMM-Newton. Highly ionized Fe 
lines, H-like transitions, and a strong underlying continuum 
in the high-resolution X-ray spectra reveal a dominant very 
hot plasma. We used three different methods to derive the 
EMD and coronal abundances and all three show a strong 
peak EM around 20-40 MK, about a dex above EMs at lower 
T (Fig.0. We compared the EMD with other giants and ac- 
tive stars (Figs.|S]and|5) to emphasize the exceptional coronal 
behavior of YY Men given its high T and EM. Such a hot 
plasma produces thermal broadening at the level detectable 
by Chandra. Indeed, we measured line broadening in several 
lines, which we in terpr eted as predominantly Doppler thermal 
broadening (Sect. ^3}. 

YY Men was about two times brighter, and possibly slightly 
hotter, during the Chandra observation than during the XMM- 
Newton observation. No evidence for flares, or small-scale 
variations down to ^ 300 s was found in the XMM-Newton 
and Chandra light curves (Sect.|5}. Nevertheless, the absence 
of variability does not imply absence of flares as the latter 
need peak luminosities > 10^' erg s"' to be detected 
with current detectors. We interpreted the shape of the EMD 
(Fig. [7} of YY Men's corona with two different formalisms. 
The first one infers the properties of coronal loops from the 



EMD shape ('Sect. l7.2.lV From the steep slope of the EMD at 
high T, we derived that most of the loops in YY Men's corona 
have their maximum T equal to or slightly above 30 MK. The 
second formalism ma kes us e of the EMD in the context of 
coronal heating (Sect. ITTT^ . We argued that a statistical en- 
semble of flares distributed in energy with a steep power law, 
dN/dE cx E-^-\ down to L^in ~ 10^^- 10^'' erg s"' could ex- 
plain the decrease of the EMD at high T and the X-ray emis- 
sion of YY Men. The steep index of the power law suggests 
that small flares could contribute most to the coronal heating 
in YY Men. 

There is a marked depletion of low-FlP elements with re- 
spect to high-FlP elements in YY Men's corona, suggest- 
ing an in verse F IP effect like in most active RS CVn bina- 
ries (Sect. m~2l i. The lack of determinations of photospheric 
abundances for individual elements except Fe does not allow 
us to determine whether the FlP-related abundance bias still 
holds when coronal abundances are compared to the stellar 
photospheric composition instead of the solar. However, a 
photospheric [Fe/H] abundance found in the literature indi- 
cates that the coronal Fe abundance is actually depleted. The 
high N abundance found in YY Men's corona is interpreted 
as a signatu re of th e CNO cycle due to dredge-up in the giant 
phase (Sect. l7.1.3> . 

The low-signal-to-noise ratios in the He-like triplets pre- 
vented us from obtaining definitive values for the electron 
densities (Sect. f74b . In addition, no significant optical depth 
was measured from line ratios (Sect. ^3}. 

In conclusion, FK Com-type giants emit strong X-rays and 
contain the hottest coronal plasmas found in the large popula- 
tion of stars with magnetic activity. Their study is important 
to understand the connection between the Sun and stars, as 
they provide the most extreme conditions (e.g., large radius, 
rapid rotation, high coronal temperature) in magnetically ac- 
tive stars. 
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Chandra and XMM-Newton Observations of YY Men 



TABLE 1 

Observation Log for Chandra and XMM-Newton 



Instrument (Mode / Filter) 


Start Stop 


Exposure (ks)" 




Chandra (ObsId 200165) 




ACIS-S/HETG 


2002 Feb 1 23:57:01 UT 2002 Feb 2 21:13:46 UT 


74.2 


XMM-Newton (Rev 334, Obsid 0137160201) 




RGSl 
RGS2 

MOS 1 (Small Window / Thick) 
MOS2 (Small Window /Thick) 
pn"^ 

OM (Fast / U) 


2001 Oct 5 12:16:31 UT 2001 Oct 6 12:50:24 UT 
2001 Oct 5 12:16:31 UT 2001 Oct 6 12:52:20 UT 
2001 Oct 5 12:23:05 UT 2001 Oct 6 12:43:09 UT 
2001 Oct 5 12:23:05 UT 2001 Oct 6 12:43:56 UT 


86.3 
84.1 
80.8 
80.8* 


2001 Oct 5 12:20:24 UT 2001 Oct 6 12:04:34 UT 


84.9'' 


"Usable exposure after filtering. 
''MOSl data only were used, 
■^pn camera off-line. 

''Fast window off-center; fixed pattern noise. 
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TABLE 2 

Spectroscopic Results from Various Methods with Chandra and XMM-Newton 







Chandra 




XMM-Newton 


Parameter 


Ml-' 


M2'' 


M3 


Ml" M2'' 



logNH(cm-2) 20.87!<;05 

C/Fe] 



N/Fe] +0-99!!! ;g 

O/Fe] +0-15!|j^ 

Ne/Fe] +0.52;g™ 



Mg/Fe] 
Al/Fe]. 



-0.03™;- 
+0.06+0" 



-2.09 

Si/Fe] +0.00+{j[^ 

S/Fe] -O-lf'tsg 

Ar/Fe] +0-42!|||l 

Ca/Fe] +"■'"'-1 42 

Fe/H] -0.52!«;i 

-Ogri(K) 6.86:0 03 

-OgTiCK) 7. 16:0 M 

.ogr3(K) 7-60!™ 

ogr4(K) ■■■ 

ogr5(K) ■■■ 

ogTeCK) ■■■ 

ogr7(K) ■■■ 

ogr8(K) ■■■ 

ogr,(K) ■■■ 

-ogrio(K) ■■■ 

-ogrH(K) ■■■ 

ogri2(K) ■■■ 

ogri3(K) ■■■ 

ogri4(K) ■■■ 

ogri5(K) ■■■ 

-ogri6(K) ■■■ 

ogrn(K) ■■■ 

ogri8(K) ■■■ 

ogri,(K) ■■■ 

.og 720 (K) 

ogTz, (K) ■■■ 

.ogEMi (cm-3) 53.95!^ ]^ 

.ogEM2 (cm-3) 5A.21^:{*. 



ogEMs (cm-3) 55.22 

EM4 . 
.ogEMs (cm"3) 

EMe (cm-3) . 
.og EM7 (cvoT^ ) . 
ogEMg (cm"3) 
.ogEMg (cm"3) ^ 
ogEMio {cwT^). 
og EMh (cm-3). 
ogEMi2 (cm"3). 
ogEMi3 (cm"3). 
og EM14 (cm-^). 
.ogEMi5 (cm"3)^ 
ogEMie (cm-3). 
og EM17 (cm-3). 
ogEMis (cm-3). 
ogEMig (cm"3). 
.ogEM2o (cm-3). 
.ogEM2i (cm"3)^ 



20-88!<!;0^ 



+0 92+"^* 
-Q-32 

+0 12+0'12 
-0.J.2 

,0 52+0 08 

i-U.JZ_„„, 

-0 03+0 11 
-0.10 
+0.O8;0-35 

-0.08!0.i6 

i-U.tO_Q 25 

+0.04+0-35 

-0 57+0:03 

6.075'' 
6.225'' 
6.375'' 
6.525'' 
6.675'' 
6.825'' 
6.975'' 
7.125'' 
7.275'' 
7.425'' 
7.575'' 
7.725'' 
7.875!" 



20.89!0.02 



20.83+^ 



5 1.06 ±0.17"= 
52.86 ±0.12=^^ 
52.01 ±0.12"^^ 
52.04 ±0.10" 
52.89±0.12'^ 
53.67±0.13'^ 
53.91 ±0.07=^^ 
53.82±0.13'^ 
53.93 ztO-H"^ 
54.46 ±0.10=^^ 
55.00±0.13'^ 
54.66 ±0.1 2=^ 
53.29 ±0.18'^^ 



+0.79!; 
+0.08!; 
+0.52: 
-0.01 



-0.09 



Statistics / d.o.f 5703/5226 5784/5219 



-0 01+010 

^•"^-0.13 
-0 09+011 

< +0.74 

+0 40+O-22 

-0 53+<' 0'' 

"■-^•^-0.07 

6.0'' 

6.1'' 

6.2'' 

6.3'' 

6.4I' 

6.5'' 

6.6" 

6.7'' 

6.8'' 

6.9'' 

7.0" 

7.1" 

7.2" 

7.3" 

7.4" 

7.5" 

7.6" 

7.7" 

7.8" 

7.9" 

8.0" 
52.07 
52.27 
52.47 
52.67 
52.83 
52.91 
52.94 
53.06 
53.23 
53.42 
53.61 
53.76 
53.91 
54.12 
54.30 
54.40 
54.45 
54.41 
54.33 
54.23 
54.08 
5947'' 



(51.67 
(52.07 
(52.45- 
(52.59 
(52.91 
(53.54 
(53.48 
(54.25- 
(54.33 
(54.08 
(53.74- 
(52.52- 
(52.92 
(53.14 
(53.10 
(53.42- 
(53.83- 
(54.15 
(54.53 
(54.47 
(54.31 



-51.87)'' 
-52.27)<i 
-52.62)'' 
-52.73)'' 
-53.04)'' 
-53.62)'' 
-53.97)'' 
-54.34)'' 
-54.23)'' 
-53.90)'' 
-52.32)'' 
-52.72)'' 
-53.07)'' 
-53.05)'' 
-53.23)'' 
-53.76)'' 
-53.93)'' 
-54.41)'' 
-54.51)'' 
-54.40)'' 
-54.16)'' 



+0.12 
+0.77 
+0.08 
+0.55: 
+0.12 



-0.11+ 



+0.18 
0.19 
_Q 45+0.29 

"■^-■-o.si 
<+0.14 

< +0.35 
-0.68!0^5 

6 91+0 03 

-flS? 

7 26+0-0' 

-0.04 
7 66+0 

/ .uu n Ai, 



53.61!*-" 
54. 16^^ 
54.8611^ 
54.68t|l 



20.82^;03 

+0.01^; 
+0.77+0-15 

+0 1 1+^^'''^ 
-AOS 
+0 51+0-10 

+"■■'1-0.09 

+0.09+0 15 



-0.13+; 



+0.14 
0.22 
_0 .54+0.30 
"■■^^-0.66 

< +0.04 

< +0.28 
_Q 64+0 04 

"■"^-0.05 

6.075'' 
6.225'' 
6.375'' 
6.525'' 
6.675'' 
6.825'' 
6.975'' 
7.125'' 
7.275'' 
7.425'' 
7.575'' 
7.725'' 
7.875'' 



52.33 ±0.34": 
53.24 ±0.25"^ 
52.52 ±0.21" 
52.56 ±0.21" 
53.16±0.21" 

53.67 ±0.24" 
53.92 ±0.23" 
54.13 ±0.24" 
54.47 ±0.22" 

54.68 ±0.21" 
54.12±0.21" 
52.85 ±0.25" 
54.11 ±0.34" 



/ 



1501/1260 1514/1260 



Note. — i) We give the EM per bin for methods 2 and 3 as de fined in Eq.|2| ii) Coronal abundances were compared to 
the solar photospheric abundances from lGrevesse &Sauval 11991 . 

"Ml and M2 obtained absolute abundances, [A/H], however we give in this Table [A/Fe] and [Fe/H] since ratios are 
more stable. The absolute abundances and their uncertainties can be calculated back, [A/H] = [A/Fe] + [Fe/H] and 
A[A/H]+ = A[A/Fe]+-A[Fe/H]-, A[A/H]- = A[A/Fe]-- A[Fe/H]+, respectively, where the + and - subscripts refer 
to the positive and negative error bars, respectively. 

''The colon indicates that the value was fixed, and corresponds to the middle of the temperature bin. 

"Obtained from the coefficients ai , . . . , av . 

'lEMD(r) for the best-fit line fluxes and its range (see text for details). As a quality measure of method 3, we give the C 
statistics after convolution of the model through the Chandra response matrices. 
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TABLE 3 

Chandra HETGS Fluxes for Lines with Excess Width. 



Ion Transition A" Flux*" 68% Conf. range'' Width'' 68% Conf. range'' Doppler width'= 
Line (A) (10-* phcm-^ s"!) (mA) (mA) 



SiXIV 6.18 18.3 (16.5-20.2) 4.4 (3.2-5.4) 2.2 

MgXII 8.42 16.0 (14.4-17.7) 5.2 (3.1-7.4) 3.2 

NeX 12.13 137 (130-143) 6.2 (5.0-7.1) 5.2 

OVIII 18.97 131 (112-153) 6.3 (2.7-9.1) 9.1 

NVII 24.78 8.0 (6.7-9.5) 4.3 (0-14.1) 12.7 



"Wavelengths from APEC 1.3.1 database. 

''Corrected tor photoelectric absorption (logiVn ~ 20.87 cm"^). 

"Obtained with AC= 1. 

''Gaussian width in excess of the instrumental profile. 
'Doppler broadening for an isothermal plasma of 40 MK. 



TABLE 4 

Chandra HETGS Fluxes Used For Density and 
Opacity-Sensitive Line Ratios. 



Ion Transition A" Flux'' 68% Conf. range" 
Line (A) (10"* phcm'^ s"') 



Si XIII 6.65 8.0 (6.7-9.5) 

Si XIII 6.69 0.9 (0-2.0) 

Si XIII 6.74 6.8 (5.5-8.2) 

MgXI 9.17 5.5 (4.1-7.0) 

MgXI 9.23 1.8 (0.6-3.1) 

MgXI 9.31 1.1 (0-2.5) 

NelX 13.45 12.9 (9.7-16.5) 

NelX 13.55 (0-1.8 

NelX 13.70 6.1 (3.4-9.2) 

OVII'' 21.60 24.4 (12.5-39.8) 

OVII'' 21.80 24.7 (12.3-40.9) 

OVn'' 22.10 5.0 (0-16.0) 

NeX 10.24 16.4 (14.3-18.6) 

NeX 12.13 137 (130-143) 

OVIII 16.00 29.2 (23.8-35.1) 

OVIII 18.97 131 (112-153) 

FeXVII 15.01 19.0 (14.9-23.6) 

FeXVII 15.26 7.2 (4.0-11.0) 

FeXVII 16.76 11.8 (7.0-17.3) 



"Wavelengths from APEC 1.3.1 database. 
''Corrected for photoelectric absorption (logA^n ~ 
20.87 cm-2). 
"Obtained with AC = 1 . 

''Very weak signal. Essentially ~ 3 counts are detected 
above the continuum for the r and i lines. 
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Fig. 1 . — Background-subtracted light curves of YY Mensae from XMM-Newton MOS 1 (a), Chandra MEGl (b), and HEGl (c) for a bin size of 500 s. 
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Fig. 2a. — (a) Chandra MEGl and HEGl spectra compared with the best-fit 3-T model (red line) found in Tab.l2l (b) Same as (a), but for 8- 14 A. (c) Same 
as {a) but for 14-26 A. The 17-26 A panels show the MEG data only. 
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Fig. 3. — XMM-Newton MOSl, RGSl, and RGS2 spectra of YY Men with the best-fit A-T model of method 1 overlaid as a thick black line. Error bars are 
plotted for MOS only, for clarity. Contributions to the (multiplied by sign[data - model]) are also plotted in the lower panel. Zoom-ins from 1 A to 4 A 
(logarithmic wavelength scale) and from 10 A to 17.5 A (linear wavelength scale) show the quality of the fits in the region where highly ionized Fe lines are 
emitted. 
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Fig. 4. — Left: Comparison between absolute abundances derived from methods 1 and 2 with Chandra (top panel) and XMM-Newton (bottom panel) data. The 
dashed line represents a 1 : 1 correlation. Absolute abundances from method 2 are slightly lower with respect to those obtained from method 1 for Chandra data 
(but not for XMM-Newton data). Right: Comparison between abundance ratios derived from method 3 and methods 1 (top) and 2 (bottom) with Chandra data. 
Excellent agreement is generally found. 
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Fig. 5. — Comparison between abundances derived from methods 1 (left panels) and 2 (right panels) with XMM-Newton and Chandra data. The dashed line 
represents a 1:1 correlation. The upper panels show the absolute abundances (relative to H), and the lower panels show the abundances relative to Fe ([A/Fe] = 
[A/H] — [Fe/H]). All abundances are relative to the solar photospheric composition iGrevesse & SauvaL1998») . 
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Fig. 6. — Abundance ratios, [A/Fe], in YY Men's corona as a function of tlie FIR For clarity, we selected abundances from method 1, and shifted the Chandra 
data points by +0.2 eV. The abundance ratios are relative to the solar photospheric composition LGrevesse & SauvaL 1998.) . The solar ratio ([A/Fe] = 0) is 
indicated by a dotted hne. 
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Fig. 7. — Emission measure distribution of YY Men obtained from the Chandra and XMM-Newton spectra for the various methods. To help comparison, we 
plot EMD(r)/Alogr, i.e., the EMDs per bin (see Table l2l were divided by the integration step (0.15 and 0.10 dex for methods 2 and 3, respectively). For the 
multi-T approach (method 1), we show the EM in each component since this method has no bin width. The y-axis, therefore, reflects this distinction (left axis for 
methods 2 and 3, right axis for method 1). The top panel shows the EMD in a linear vertical scale to emphasize the dominant very hot plasma, whereas the lower 
panel uses a logarithmic vertical scale to reveal the weak, but required plasma at lower temperatures. 
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Fig. 8. — Comparison of the HMDs obtained in tliis paper for YY Men (method 3, using the definition in Eq.|2j and those derived for the giants Capella, 
31 Com tArgiroffi et al_2003.:.Scelsi et al_20Q4.) . Note the much flatter HMD, the higher temperatures, and activity level of YY Men. 
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Fig. 9. — Comparison of the results of 4-7 fits of YY Men XMM-Newton spectra with those obtained by 'Audard et al.'f2003") for a sample of RS CVn stars of 
different activity levels. Also s how n here is a 3-7 fit o f the sin gle rapidly rotating giant 31 Com ( Scelsi et al. 2004). Fits for two other FK Com-type stars (V1794 
Cyg and FK Com) from lSondoiri jSOPjl) and lGondoin et alj 120021) are also shown. The lines connecting the data points of each star are an aid to the eye and 
have no physical meaning. 
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Fig. 10. — Radiative cooling curve in units of erg cm' s ' for T = 1 — 100 MK, based on the APEC 1.3.1 database (for photon energies from 0.01 keV to 50 keV, 
i.e., 0.25 A to 1240 A) and for best-fit XMM-Newton abundances (Ml; Tab.l2l. The total radiative loss function (solid) is the sum of the continuum contribution 
(dash-dotted) and of the Une contribution (dashed). For comparison the total radiative loss for solar abundances iGrevesse & Sauva]»1998.) is shown as a dotted 
line. 
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Fig . 11. — Radiative loss time, Tr(T) = 3kT / (/Je A(r)), in units of days derived from the radiative cooling curve, A(r), in YY Men's corona as shown in Fig. 
Two different plasma electron densities, itc = l(f and lO'" cm"^^, were used. As mentioned in Sect. 17. 4j there is no evidence for higher densities in YY Men. 
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Fig. 12. — (Top panel): Line photon luminosity distribution for Ne X Lya (solid) and O VIII Lya (dotted). (Bottom panel): Weighted average width 
distribution for the same lines (see text for details). 
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Fig. 13. — (Left) Confidence map for the Ne X Lya fine. Confidence contours with AC = 1 and AC = 2.71 are shown in red (inner dotted contour) and green 
(outer sohd contour), respectively. The loci of the best-fit parameters for 10,000 Monte-Carlo simulations of an emission line with an instrumental profile are 
shown as cyan-colored dots. (Right) HEG line profile of Ne X Lya in the negative (top) and positive (bottom) order spectra. The best-fit Gaussian profile is 
shown as a red line (wide solid line), whereas the best-fit with the instrumental profile is shown as a green fine (narrow dashed fine). Residuals from the fit with 
a Gaussian profile are shown in sub-panels. 
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Fig. 14.— Similai- to FigurefHlbut for Si XIV Lya. 



